%根据构型config计算体系能量
function en = en_cal(config)
    N = length(config);
    en = 0;
    for i=1:(N-1)
        for j=(i+1):N
            en = en + LJ_pot(config{i},config{j});
        end
    end

    return
end

%取原子直径为0.88sigma，截断半径为3.0sigma
%pt1,pt2均为1*3的数组，记录原子中心在三维直角坐标系下的3个坐标分量
function pot = LJ_pot(pt1, pt2)
    %根据pt1,pt2两点的坐标计算两点距离
    rvector = pt1 - pt2;
    r_square = sum(rvector.^2);
    r = sqrt(r_square);
    pot = r^(-12)-r^(-6); 
    
%     if(r < 0.88)
%        pot = 0.88^(-12)-0.88^(-6);
%     elseif(r > 3)
%         pot = 0;
%     else
%         pot = r^(-12)-r^(-6);   
%     end

    return
end